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Abstract 

Renyi entropies S g are useful measures of quantum entanglement; they can be calculated 
from traces of the reduced density matrix raised to power q, with q > 0. For (d + 1)- 
dimensional conformal field theories, the Renyi entropies across S d ~ l may be extracted from 
the thermal partition functions of these theories on either (d + l)-dimensional de Sitter space 
or K. x M d , where M. d is the rf-dimensional hyperbolic space. These thermal partition functions 
can in turn be expressed as path integrals on branched coverings of the (d + l)-dimensional 
sphere and S 1 x M d , respectively. We calculate the Renyi entropies of free massless scalars 
and fermions in d — 2, and show how using zeta-function regularization one finds agreement 
between the calculations on the branched coverings of S* 3 and on S* 1 x M 2 . Analogous 
calculations for massive free fields provide monotonic interpolating functions between the 
Renyi entropies at the Gaussian and the trivial fixed points. Finally, we discuss similar 
Renyi entropy calculations in d > 2. 
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1 Introduction 



The Renyi entropies [TJ|2) have recently emerged as powerful diagnostics of long-range en- 
tanglement in many-body quantum ground states in d > 2 spatial dimensions [3 15 . The 
Renyi entropy S q is defined as 



1-q 



log tr p q , q > 



where p is the reduced density matrix obtained after tracing over the degrees of freedom 
in the complement of the entangling region. Numerical evaluations of the entropies have 
allowed characterization of many distinct types of ground states: gapped states with topo- 



logical order 4148 , quantum critical points 9 , Fermi liquids 10,11 , non- Fermi liquids 12 



and Goldstone phases with broken symmetries 13-15 . Indeed, it appears that the Renyi 



entropies carry distinct signatures of the known quantum many-body states and are also 
amenable to evaluation by convenient algorithms. 

Nevertheless, there are few known analytic results for Renyi entropies with general q. In 
d — 1, corresponding toD — d+1 — 2 space-time dimensions, for gap less conformal states, 



the Renyi entropies of an interval of length R are given by 16 18 



S a 



\og(R/e) , d=l 



:i-2) 



where e is a short-distance cutoff, and c is the central charge of the conformal field theory 
(CFT). For odd d > 1, in studying Renyi entropies across the (d — l)-dimensional sphere of 
radius R, one again finds terms logarithmic in R. Their coefficients as a function of q were 



calculated in 19 21 for a massless real scalar field. In the case d = 3, the result is 

(l + g)(l + <^ 



s. 



scalar 



a 



360g 3 



logOR/e) 



;i.3) 



where a is a non-universal constant dependent upon ultraviolet details, and e is again a 
short-distance cutoff. For q = 1 this result agrees with a direct numerical calculation of the 



entanglement entropy 22 . We note for completeness that logarithmic terms also appear in 

- l/q)R d ~ 1 InR, but these modify 



10 



12 



Sq ~ 



the Renyi entropy of a d > 1 Fermi liquid 
the leading 'area law' term and so are much stronger than those in (1.3). 

In d — 2, there is a simple known result for gapped topological states, such as those 
described by an effective U(l) Chern-Simons gauge theory at level k. The Renyi entropies 



1 



are independent of q and given by [7[[9| [23[|25] 

R 1 



S q = a---hxk, (1.4) 
e 2 



where R, a, e are as in (1.3). Note that there is no logarithmic dependence on the size R in 
the universal term. We reproduce this result in Appendix [Cj 

Generally the structure of Renyi entropies in conformal field theories (CFTs) in d = 2 



is more complicated than the structure in d — 1 26 -28 . In d = 2 the leading term ~ R/e 
is UV divergent and not universal. However, the subleading /^-independent term is finite 
and universal (in even d such universal terms have no logarithmic dependence on R). Unlike 
in d — 1, its dependence on q is not simply through a general g-dependent normalization 
factor. More generally, we can deviate away from the CFT with a relevant operator which 



generates a mass scale m; then the Renyi entropies obey 27 



R 

S q = a- - S q (mR) , d=2, (1.5) 



where S q (x) is a universal function, which in turn obeys 



r q x as x — > oo . 



S q (x) = { (1.6) 
7 g as x — y , 

where r q and 7 g are universal numbers. The number r q describes the shift in the linear R 
dependence of S q due to the presence of the mass scale m, while 7 g measures the universal 



contribution of the CFT to the Renyi entropy, similar to eq. (1.4). Ref. 27 obtained results 
for an infinite cylinder of circumference R divided along a circular boundary: explicit results 
were obtained for all q for both j q and r q for the CFT of a free scalar field, while only r q 
was determined for general q in the large N limit of the Wilson-Fisher CFT of ^-component 
interacting scalars. 

All our subsequent discussion of Renyi entropies will be restricted to conformal field the- 
ories in which m — 0. Also, we will not display the leading non-universal term proportional 
to R, and so we write simply 

S q = - lq . (1.7) 



In the limit q — )■ 1, the Renyi entropy (1.1) becomes the more familiar quantum entan- 



glement entropy (see (17||29||30| for reviews and references to earlier work). While in general 
its calculation is quite difficult, simplifications do occur for conformal field theories with 



2 



some particular geometries separating the entangling regions A and B. For example, in the 
case of entanglement across a sphere S^ -1 in d flat spatial dimensions, the use of conformal 
mappings reduces the calculation of entanglement entropy to that of the thermal entropy 



evaluated when the spatial geometry is taken to be the hyperbolic space M d 19 ,28 31 
For the purpose of calculating the entanglement entropy, the temperature is taken to be 
To = l/(2irR), where R is the radius of curvature of M. d . For this special temperature, the 
required Euclidean path integral calculation on S 1 x M. d is related by a further Weyl trans- 
formation to that on S d+1 |31|. For example, in d = 2 the entanglement entropy between a 
disk and its complement is given by — F = log \ Z\, where Z is the Euclidean path integral 
on S 3 . The latter quantity can be calculated in any 3-dimensional field theory with J\f > 2 



supersymmetry using the method of localization 32 35 . It has also been calculated in some 
simple non-supersymmetric CFTs, such as free field theories [l9|[28|[36|[37] , the Wilson-Fisher 
fixed point of the O(N) model for large N [37] , and conformal gauge field theories with large 
numbers of flavors |38|. 



A separate approach to the computation of entanglement entropies [39] (see also 40 ,41 



relies on the AdS/CFT correspondence 42 -44 . This approach has been used to calculate 
entanglement entropy across S' d ~ 1 in CFTs that possess a dual description with weakly curved 



gravitational backgrounds 31 45 . All these results have led to conjectures that the 3-sphere 
free energy F, or equivalently minus the entanglement entropy across a circle, decreases 
along any RG flow and is stationary at the fixed points [3l|[37|[45| [46 1 f] 

The conformal methods for calculating the entanglement entropy S\ have been extended 
recently Jl9|[48|[49| to allow computation of Renyi entropies S q in CFTs for spherical entan- 
gling geometries [l9j|49]. It turns out that now one needs the thermal free energy of the 
theory on M. d evaluated at temperatures T Q /q. It is useful to define 



Fq = "log \Z, 



91 ' 



where Z q is the Euclidean path integral on S 1 x H , with S l of circumference 2irRq. This 
geometry can be further conformally mapped to a certain g-fold branched covering of S d+1 
analogous to the one described in the next section for d = 2. While this branched covering 
has curvature singularities, we will show that the path integral can often be simply calculated 
on this background, and the results agree with those on S 1 x M d . The Renyi entropies can 



1 After the original version of this paper appeared, a proof of the F-theorem was presented in 47 . 
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then be simply expressed as 19 49 



s q = q:F ; Tq . (1.9) 



In 49 these methods were applied on the gravity side of the AdS/CFT duality, where the 



free energy is read off from that of a certain topological black hole geometry. The goal of 



our paper is instead to continue along the lines of 19 and consider applications of these 
methods to free field theories, such as the conformally coupled scalar field and the massless 
Dirac fermion in d = 2 (D = 3). 

We will present calculations both on S l x M 2 and on the g-fold branched coverings of S 3 
and demonstrate their consistency. For simplicity we only consider integer g-fold branched 
coverings of S 3 (some further details may be found in Appendix A), though in principle the 
partition function on that space can be computed for any q > 0. We consider the conformal 
scalar field in section 3, and the massless Dirac fermion in section 4. In section 5 we show 
that analogous calculations for massive free fields provide monotonic interpolating functions 
between the Renyi entropies at the Gaussian and the trivial fixed points. In the Discussion 
section we comment on extensions to d > 2. 



2 Methods for computing the Renyi entropy in CFT 



In 31 ,49] it was explained how one can use conformal symmetry to express the Renyi entropy 
(1.1) in terms of partition functions of the CFT on certain curved manifolds. The curved 
manifolds relevant for d = 2 are multiple branched coverings of the three-sphere and S 1 xM 2 , 
where H 2 is the two-dimensional hyperbolic space. One can see how these spaces arise by 



following the arguments presented in 31 49 . With the help of a Weyl transformation one 
can map the reduced density matrix on a disk of radius R in IR 2,1 to 



e ~2nRHt 



(2.1) 



where Ht is the Hamiltonian generating time translations in either: (A) the static patch of 
three-dimensional de Sitter space of radius R with the metric 

ds\ = - cos 2 6dt 2 + R 2 [d0 2 + sin 2 6d<f) 2 ] ; (2.2) 



4 



or (B) R x H 2 , with the metric 



ds 2 B = -dt 2 + R 2 [dr] 2 + sinh 2 rj dcfr 2 ] . 



(2.3) 



That the two spaces above differ just by a conformal rescaling can be seen by writing ds\ 
cos 2 9 ds 2 B , with the identification sinh rj = tan 9. 



Using the definition of the Renyi entropies in eq. (1.1) and the density matrix (2.1), one 
finds 



1 , tie -^RqH t 
S n = -, log 



qT x -T q 2uRq (Ft - F q 



where we made the definitions 



1-q 



l-q 



(2.4) 



Z q = tr e~ 27rRqnt =e~ T \ 



T q = 2nRqF t 



(I ! 



(2.5) 



which hold for both of the mappings (A) and (B) introduced above. 

The quantity Z q can be interpreted as the partition function at temperature T = l/(2irRq) 



in either de Sitter space in the case (2.2) or R x H in the case (2.3), and F q = F(T) is the 



corresponding thermal free energy. As is standard in thermal field theory, Z q can be com- 
puted as a Euclidean path integral after the Euclidean time direction has been compactified 
into a circle of length (3 = 1/T = 2nRq. When we use the mapping (A) to de Sitter space, 
the Euclidean metric is 



ds 2 A = R 2 [cos 2 Bdr 2 + d6 2 + sin 2 6d(j> 2 ] 



(2.6) 



where the ranges of the coordinates are < r < 2ixq, < < 2tt, and < 6 < tt/2. 
Similarly, when we use the mapping (B) to R x HI 2 , the Euclidean metric is 



ds 2 B = R 2 [dr 2 + dr] 2 + sinh 2 r] d(j) 2 ] , 



(2.7) 



where the ranges of the coordinates are < r < 27T q, < cf) < 2tt, and < rj < oo. 



Like their Minkowski counterparts (2.2) and (2.3), these Euclidean spaces also differ from 
each other just by a conformal rescaling, as can be seen from writing ds\ = cos 2 9 ds 2 B with 
sinh 7/ = cos#. 



Before we examine the spaces (2.6) and (2.7) in more detail, let us comment on the case 



5 



q = 1, where eq. (2.4) seems poorly defined, and rewrite this equation in a different way. 



Taking the limit q — > f in (2.4), one obtains 



Si 



dF(T) 



dT 



-S therm (l/(27ri?)) 



T=l/{2nR) 



(2.8) 



Stherm{T) being the thermal free energy at temperature T. In general, the definition of the 
free energy is 



F(T) = E(T) — TS thcim (T) 



(2.9) 



where E(T) = tr(pq'Ht) is the total energy. Precisely at T = l/(2irR), or equivalently at 
5 = 1, the manifolds (2.6) and (2.7) are conformally equivalent to M 3 , which implies that for 



either of these manifolds we must have E(1/(2ttR)) = 0. It follows immediately from (2.8) 
and (2.9) that Si = —2-kRFi = —J-\. A simple calculation then shows 



S q — Si 



•Fq — •/ 1 

9-1 



2irR 



q Fq - F\ 
q-1 



(2.10) 



Let us now comment on the features of the spaces (2.6) and (2.7). The space in eq. (2.6) 



is not smooth unless q = 1, in which case it reduces to the round three-sphere of radius R. 



One way to see this fact is to think of S 3 as the sphere in C 
We can write 



2 ~ 



;iven by \z±\ +\z%\ 



R 2 



Zi = Rcos9e l 



Z2 = R sin 9e 



icj> 



(2.11) 



where 9 ranges from to n/2, and r and <p range from to 2n. These coordinates are the 
same as those appearing in (2.6), as one can easily check that when q = 1 the standard 



line element in C takes the form (2.6) under the mapping given by (2.11). Making r range 



from to 2nq one then obtains (2.6) for arbitrary q. The description of (2.6) in terms of 



the complex coordinates Z\ and Z2 is particularly useful when q = 1/p for some positive 
integer p. In this case the space (2.6) should be identified with the orbifold S 3 /Z p , the Z p 

? 2nl /p. On the other hand, when q takes values in the 



equivalence being given by z\ 



Z\e 



positive integers, (2.6) describes a g-fold covering of S 3 (henceforth denoted C q ) branched 
along the circle located at 9 = 7r/2. 

One way of defining the spaces C q , and really the way we will be thinking about them 
in this paper, is in terms of the normalizable functions (or more generally sections of spinor 



6 



and vector bundles) that can be denned on them. For instance, when q = 1 a basis for the 
Hilbert space of normalizable functions consists of polynomials in z\, z%, and their complex 
conjugates; when q = 2, we are also allowed to have Jz\ or t/zf times polynomials in z\ 
and Z2-, but not yfz^ and y/z% times such polynomials, etc. As we will see, the set of allowed 
functions and sections of spinor and vector bundles on C q is one of the central ingredients of 
our computations. 



One of the advantages of working on the spaces (2.6) as opposed to (2.7) is that these 



spaces are compact, and after removing the UV divergences (for example by zeta- function 
regularization) the free energies F q are finite. On the other hand, the hyperbolic cylinder 



(2.7) (henceforth denoted H q ) is non-compact, and regularization of UV divergencies yields 
a finite free energy density. The free energy itself is proportional to the volume of H 2 , which 
is infinite and requires further regularization. The proper regularization of this volume [50] 
uses a hard cutoff at some value rj = r] : 



no 



Vol(ff) = 2tt / dr] sinhr? = 2tt 



o'/0 



-1 + 



2 2 

Taking the finite part of this expression, we should identify the regularized volume as 



(2.12) 



Voi(e 2 



-2vr. 



(2.13) 



3 Renyi entropies for free conformal scalars 

Among the simplest conformally invariant field theories in d + 1 = D > 2 are free scalar 
fields cf) conformally coupled to curvature. The action for such a complex field on a manifold 
M is 



S 



d D r^/g(r) 



\0„ 



+ 



D 



A(D - 1) 



n |0| 



(3.1; 



where 1Z is the Ricci scalar of M. In this section we will set D = 3 and work on the spaces 
C q and H q introduced in the previous section. These spaces have constant scalar curvature 
7Z, so the path integral yields 



J= = tr log -A + 



K 



(3.2) 
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where A is the Laplace operator. We will first calculate T when M = C q (TZ = 6) and then 
reproduce the same results by taking M — H q (1Z — —2) (from here on we set the radius 
R = 1). 



In general, 1Z = —d(d — 1) on S 1 x M d . Thus, the conformally coupled scalar (3.1) has 



m 2 = — n = - {d 1)2 . (3.3) 

Ad 4 

This means that the constant mode on S 1 saturates the Breitenlohner-Freedman (BF) bound 



51 ontf. 



3.1 Conformal scalars on the g-fold branched covering of S 3 

The computation of J 7 on C q starts with the diagonalization of the Laplacian on this space 
(in this section we take q to be a positive integer). We will show that the Laplacian has 
eigenvalues 

A n = -n(n + 2) (3.4) 

for every n E N/g0with de generacy 

f(n + l) 2 , ifnGN, 
9n = < (3.5) 
[(ib + l) (A; + 2) , ifn = fc + 2, fc,pGN, l<p<q. 

This result can be proven by examining the differential equation satisfied by an eigenfunc- 
tion of the Laplacian. The space C q has two U(l) isometries corresponding to shifts in r and 
<ft, so the eigenfunctions of the Laplacian can be assumed to be of the form / '(#)e" iv ' r+rm <^ 
for some m T , m^, and f(9). The eigenfunction equation for / with eigenvalue A is 

izre + *k) m = xt ®- (3 - 6) 

The solution regular at 9 = n/2 is 
f{9) = (cos0) |mT| (sin0) m ^ 

(3.7) 



x 



/ 1 + \m T \ + - Vl - A 1 + \m T \ + + y/1 - A . 2 A 

2-Pi ( 2 ' 2 ' T ' C ° S / ' 



2 We use N to denote the set of non- negative integers (including zero). 



S 



Regularity at 9 = implies 



A = — n{n + 2) , n = \m T \ + \m^\ + 2a , aGN. 



(3.8) 



Note that n itself is not required to be an integer. As discussed towards the end of section [2j 
the allowed values of m r and are related to the periodicities of r and 4>, respectively: the 
period of <j) is 2ir and that requires G Z; the period of r is 27rg, and the condition for 
having well-defined eigenfunctions is m T G Z/g. It follows that n G N/g, as claimed above. 
It is then straightforward to show by induction that the number of ways one can choose m T , 



m^, and a so that n = \m T \ + |m^| + 2a is given by (3.5). 



According to (3.2) the free energy is 



k=l 



k 2 log 



<J — 1 oo 

y;y>(*+i)]o 6 



p=l k=l 



k + 



V 



(3.9) 



This expression is divergent and requires regularization. The first step in obtaining a reg- 

k + 2 ' " 



ularized expression is to write log 



log 



+ log 



* + !-§ 



after a rearrangement of the terms in the sums one obtains 



and 



J^ = ^(2A; 2 -2A; + l)loj 



k=l 



q—1 oo 

p=l fe=l 



P 



(3.10) 



We can now use zeta-function regularization to write 



E 

k=l 



2k 2 -2k + l 



A- 2 



g— 1 oo 

P-i *.i (* + f - 3 



s=0 



(3.11) 



The sums appearing in (3.11) are easily evaluated in terms of Hurwitz zeta-functions: 



•Fg = -2 



C'(-2, 1/2) + ^'(0,1/2) 



p=l 



p 1 
q 2 



q 2) V q 2 



(3.12) 



where the derivative acts on the first argument of the Hurwitz zeta-function. 



It is also useful to write T„ in a more elementary form using the identities (B.2) presented 



9 



in Appendix [Bj With the help of these identities one can show 



q ~ (-1)" 1 ~ cot(n7r/g)(-ir 1 y> a n 

q 2tt 2 ^ n 3 2tt^ n 2 2q ^ n ' (3-13) 

n=l n=l n=l 

q\n q\n 



where 



•1)™ csc 2 (mr/q) ifgfn, 
•1)"2!±2 ifgjn. 



(3.14) 



Note that all the sums in eq. ( |3.13[ ) are convergent: the first two are absolutely convergent 

l+2)q 
-mq+\ 



and the third converges because Yur^mg+i a n = for any m E N. 



The first few free energies are 
log 2 3C(3) 



4 8tt 2 
log 2 , C(3) 



4 87T 2 ' 

T _ log2 C(3) V>i(l/6) + ^i(l/3) - ^(2/3) - ^(5/6) (3 ' 15) 

3 4 24tt 2 + 72V3n 
_ 51og2 C(3) g 

4 16 32vr 2 + 2tt ' 

where ipi(x) is the digamma function, and G is the Catalan constant. It is also not hard 
to find the large q asytmptotics of J- q . When q is large, we can approximate cot(mr/q) w 



csc(mr/q) ~ q/(nn) in (3.13), and so 



ra -^E^? = ^ at lar § e « ■ ( 3 - 16 ) 



n=l 



3.2 Conformal scalar on the hyperbolic cylinder 



In this section we calculate T q on the hyperbolic cylinder H q introduced in eq. ( 2.7[ ) (now we 



do not restrict q to be an integer). In anticipation of section [5] below, we start by considering 
a massive complex scalar field on H q with the action 

S = [ d*r^W) [\d,<P? + M 2 |0| 2 ] . (3.17) 
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The conformally coupled scalar (3.1) corresponds to M 2 = —1/4, and we will shortly spe- 



1S 



cialize to this value of M 2 . It is convenient to define m 2 = M 2 + j, so that the confor- 
mally coupled scalar has m 2 = 0. Since on the M 2 of unit radius the BF bound 51 
M 2 > M 2 BF = -i, we note that m 2 = M 2 
The free energy on H q is 



M 2 BF . 



JF q {m) = trlog 



-A \-m' 

4 



(3.18) 



The operator under the log is diagonalized by wavefunctions of the form f(j], (j))e inTT ' q for 
integer n T and for / being an eigenfunction of the Laplacian on M 2 . Since M 2 is not compact, 
the spectrum of the Laplacian on it is continuous, consisting of eigenvalues A + 1/4 for A > 
with the density of states 



V(X)d\ = Vq1 ^ M ^ tanh(7rv / A)rfA . 

47T 



(3.19) 



(see for example 52 -54 ) . The eigenvalues of — A + m 2 — 1/4 on H q form towers of continuous 
states for each n T , each tower consisting of eigenvalues 



A + -£ + m 2 



with the same density of states as (3.19). The free energy (3.18) becomes 



JF q {m) 



E 



n T =— oo 



nz 



dXV(X) log A + -f + m 



We can use (B.3) to compute the regularized sum over n and we obtain 

roc 

jr g ( m ) = / dXV{X) [2 log ^1 - e - 2nqVX ^J + 27cqVX + m 2 



(3.20) 



(3.21) 



(3.22) 



We will use this formula for arbitrary m > in section [5j but for now we set m = 0. 



The integral in eq. (3.22) is divergent, but the divergence is in the temperature-independent 



vacuum energy given by the second term in (3.22). This divergence can be regulated with 
zeta-function regularization or equivalently by subtracting the 1R 3 free energy density 



q- 



Voi(e 2 



dX\f\~. 



(3.23) 
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Performing the integral of the difference between the second term in (3.22) and (3.23), and 
using Vol(H 2 ) = — 27r as explained at the end of section [2J we obtain 



/oo , . 

dAtanh(TrVX) log (l - e - 27nzVX J + q 



303) 
4vr 2 



(3.24) 



This expression is convergent for any q > and it represents one of our main results. 



While we are not aware of a general formula for the integral in (3.24), when q is an 



integer the expression (3.24) evaluates to the free energies (3.12) — (3.15) we found in the 



previous subsection from the g-fold branched covering of S 3 . For example, when q = 1, one 
can evaluate the first integral: 



/oo 
rfAtanh(vrv / A) log (l - e' 2w ^J 



log 2 9C(3) 



(3.25) 



which, when combined with (3.24) yields precise agreement with (3.15): 

log 2 3C(3) 



8vr 2 ' 



(3.26) 



Similarly, 



ft- ft 



d A t anh(,VA)log 1+ta " h( ^ ) + 3C(3) - C(3) 



2tt 



2 - 



(3.27) 



also in agreement with (3.15), etc. In the limit of large q, the second term in (3.24) clearly 



dominates, and its value agrees with (3.16) 



As another consistency check, we can demonstrate by explicit computation that the total 
energy at temperature Tq = 1/(2tt) vanishes. Indeed, in general 



E = F -T— = - 1 

dT 2n dq 



(3.28) 



Differentiating our general formula (3.24) with respect to q, we obtain 



E(T ) 



1 &T„ 



2n dq 



d\ 



I + e 27rv / A 87T 3 



(3.29) 



confirming the conceptual argument used in section [2] that yielded the same result. 
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3.3 Results for the Renyi entropies of a complex conformal scalar 



Now that we have checked the agreement between the calculations of T q — — log \Z q \ on C q 



and H q , we can use (1.9) to calculate the Renyi entropies. Some of the results for integer q 
are listed below: 



Sl = = _!!fH + m « -0.1276 



S2 — J~2 ~ 2J 7 1 
5.3 



8tt 2 
log 2 | 7C(3) 



4 



-0.0667 



1 t: _ 3 + 13 C( 3 ) , ^i(l/6) + ^i(l/3) -^(2/3) -^(5/6) ^ 

2 3 2 1 4 24vr 2 + UAV^ti ~ ' 



5 4 - ^4 - 3^1 



11 log 2 49C(3) g 
48 96tt 2 6tt 



S a 



.jFj + lim 5? 

g^oo q 



4 



-0.0481 



log 2 9C(3) 

+ -0.0363. 



(3.30) 



At large g we find that 



S g = S 00 (l + ^ + 0{l/^)J . 



(3.31) 



A plot of S'q is presented in Fig. [Tj We note that, as in the holographic calculations of [49], 
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-0.08 
-0.10 
-0.12 



Figure 1: The Renyi entropies S q for the complex conformal scalar field. Note that the 
function S q is a monotonic function of q. The black dashed line is the asymptotic value S^. 



S q is a monotonically increasing function of q. It would be interesting to compare (3.30) 



with direct numerical calculations similar to those performed for d = 3 in 22 
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4 Renyi entropies for free fermions 



Another simple conformal field theory is that of free fermions. On a manifold M of dimension 
D, the action is 

S = J d D r^{r)i,{iIp)ijj, (4.1) 

where ilp is the Dirac operator. The factor of i appears in the action such that the operator 
ilp is Hermitian, and hence it has real spectrum. The free energy one obtains from integrating 
out ip is 



7 = -trlog(0) . 



(4.2) 



We will now compute J 7 in the case M = C q and then in the case M 



H q , in D 



3. We 



then use eq. (2.4) to extract the Renyi entropies. We will only study Dirac fermions, which 
are two-component complex spinors in D = 3. 



4.1 Free fermions on the g-fold branched covering of S 3 

In order to compute J 7 on C q , one first has to diagonalize the Dirac operator. The diag- 
onalization can be done by examining the eigenvalue problem in differential form, as was 
done for the Laplace operator in section [3} The analysis in the case of the Dirac operator is 
more involved that that in section [3j because one now has to deal with spinors as opposed 
to functions, and we leave the details for Appendix A. 2 Here let us just quote the result: 



the eigenvalues of the Dirac operator on the qth branched covering C q of S 3 form q towers 
parameterized by an integer p with < p < q. For each p the eigenvalues are 

V 1 
q 2q 

with degeneracy 



±[k + l + - + ^-) , keN, (4.3) 



(h + l)(k + 2) (4.4) 

for each choice of sign. For example, when q = 1 we have only one tower of eigenvalues 
k + 3/2; when q = 2 we have two towers, one consisting of/c + 3/2 — 1/4 and one of 
A; + 3/2 + l/4, etc. The various towers are symmetric about A; + 3/2, but the numbers fc + 3/2 
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are part of the spectrum only for odd q. 



The free energy following from the eigenvalues (4.3) and degeneracies (4.4) is 

q— 1 oo 



y oo / -, \ 

^ = - 2 EE^ + i)iog ( k + - q + y q ) 



(4.5) 



p=0 k=0 



This expression is of course divergent, but after zeta-function regularization it can be put in 
the form 



9-1 
p=0 



1 q 2q) \ 2 q 2qJ ^ V ' g 2g 



1 p 1 



(4.6) 



where the derivatives act on the first argument of the Hurwitz zeta-function as before. While 
this is our final answer for the free energy of fermions on C q , it may be instructive to put 



it in a form similar to eq. (3.13) from the boson case by using the zeta-function identities 



(B.2): 



~ 3£(3) 1 csc(n7r/g) 1 co 



i >ir 2 q 2 2tx * — ' n z 

71=1 

q\n 



cos(nir/q) 2q 2 + 1 

+ ~ o log 2 



2q nsin (mr/q) 12q 

qjn 



(4.7) 



These sums can be evaluated in terms of polygamma functions, but the resulting expression 
is not very enlightening. 

The first few T q can be obtained from evaluating the sums: 



~ 3C(3) log2 
J~\ — r, n H ; — 



~ _3C(3) 31og2 G 
2 32vr 2 16 2tt ' 
~ C(3) , Vi(V3) , log 2 2tt 



(4i 



24tt 2 



V3 ' 



As in the case of bosons, it is also possible to find the large q asymptotics of T q by making 
a small angle approximation inside the trigonometric functions: 



1 n 2 t — i n o 

n=l 



oo 1 

Y- 



# 9 atl^ge,. 



7T Z 



(4.9) 
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4.2 Fermions on the hyperbolic cylinder 

Like for bosons, we can reproduce the results of the previous subsection from a computation 
on the space H q . Let us start with the action for a massive Dirac fermion of mass m and set 
m = later on. The action is 



S = d 3 r a/ g(r) ip (iJp + im) ip . 



(4.10) 



The free energy obtained by integrating out ip is 



T q (m) = — tr log [ilj) + im) . 



(4.11) 



In order to calculate J- q {m) we need to diagonalize the operator ilp + im. 

The space H q has an isometry corresponding to shifts in r that commutes with the 
operator + im, and this isometry allows us to assume that the eigenfunctions of this 
operator are of the form e mT l q times a spinor on H 2 . It is not hard to show that for each n 
there is a continuous spectrum of eigenfunctions parameterized by a positive real number A, 
with eigenvalues 



Tl 

±\l X 2 + ^ + im . 



(4.12) 



The density of states for the continuous parameter A > is [53 , 54 



(4.13) 



7T 



We require the fermions to be antiperiodic on S 1 , so n is allowed to take only half odd- integer 
values. The free energy is then 



JF q {m 



dXT)(X)\og 



A 2 + % + m 2 



(4.14) 



We can use (B.3) to do the sum over n: 



F a (m) = -- 



d\V(X) 



2 log ( 1 + e" 2f9v ^ ) + 2tt(7\/A 2 + m 2 



(4.15) 
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We now set m = 0. With Vol(H 2 ) = — 2ir, we have 

POO 

T q = f q (0) = / dX Acoth(TrA) [2 log (l + e~ 2nqX ) + 2vrgA] . (4.16) 
Jo 

The first term leads to a convergent integral, while the second term needs to be made finite 
by either zeta-function regularization or by subtracting the flat space free energy density. 
Both of these procedures yield the same answer, so we will use the latter: 



2nq / dX A 2 coth(vrA) = 2vrg / dX A 2 (coth(vrA) - 1) = q 
Jo Jo 



TV 



2 ' 



(4.17) 



We can now write down the free energy of a massless Dirac fermion in terms of a conver- 
gent integral: 



POD 

7 q = 2 / dXX coth(TrA) log (l + e~ 27r9A ) + q 
Jo 



C(3) 



7T 



2 " 



(4.18) 



Along with the corresponding result for a conformally coupled scalar in eq. (3.24), eq. (4.18) 
is another main result of this paper. 



We can check in particular examples that eq. (4.18) reduces to eq. (4.7) that was derived 
on the g-fold branched covering of S* 3 . For example, when q = 1 we have 



T x = 2 / dXX coth(TrA) log (l + e~ 27rA ) + 



-2^ , C(3) log 2 , 3C(3) 



+ 



8tt 2 



(4.19) 



which is in agreement with eq. (4.8) and 37 . Also, when q is large the first term in (4.18) 



becomes negligible compared to the second term, in agreement with the large q asymptotics 



provided in (4.9). 



As in the case of bosons, we can check explicitly that the total energy at temperature 
Tq = l/(27r) vanishes: 



E(T ) 



1 &T„ 



2n dq 



9=1 



1 + e 



2ir\ 



7T- 



(4.20) 



This result confirms the argument based by conformal mapping to 1R 3 presented at the end 
of section [2l 
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4.3 Results for Renyi entropies of a Dirac fermion 



Using (2.10) we can now calculate the Renyi entropy of a free massless fermion. Here are 



some particular cases: 



Si = -Ti 



log 2 3C(3) 



-0.21896 



57T 



16 32tt 2 2tt 



-0.15076 



63 " Y z ~ Y 1 



log 2 13C(3) + Vi(l/3) 



S 

At large q 



— T\ + lim 



4 24vr 2 6\/37r 9^ 
f q _ log 2 , 5C(3) 



8tt 2 



-0.09717. 



S q = S oc (l + - + 0(l/q 2 



-0.13157 



(4.21) 



(4.22) 



A plot of S q is presented in Fig. [2j We note that, like in the scalar case, S q is a monotonic 
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-0.22 
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Figure 2: The Renyi entropies S q for the massless Dirac field. Note that the function S q is 
a monotonic function of q. The black dashed line is the asymptotic value Soo- 

function of q. 



5 Massive free fields 



In this section we study more closely the theories of massive free scalar and Dirac fields on 
H q . If the mass of such a theory is denoted by m, where m = corresponds to the conformal 
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case, then we calculate the functions 



l-q 



(5.1) 



where J r q {m) is the free energy of the theory on H q , with radius set to unity 

While by construction f g (0) = S q , for generic m > the function f q (m) is not equal to 
the entanglement entropy across a two-sphere in three flat dimensions. The mapping of the 
Renyi entropy in (^-dimensional flat space across an S"^ -1 entangling surface to a calculation 
on H q is only valid at conformal fixed points. We will see, however, that in the limit m — >■ oo 
each f q (m) falls off to zero exponentially, as is expected of the long range entropy of a massive 
theory. The function f q (m) can therefore be thought of as an interpolating function for the 
qth Renyi entropy, between the Gaussian and the trivial conformal fixed points. We observe 
for free scalar and Dirac fields that this function is monotonic along the RG flow. 

5.1 Free massive scalars 



The action for a complex massive free scalar field on H q is given in equation (3.17). The 



free energy for this theory on H q can be found in equation (3.22). The second term in 



equation (3.22), corresponding to the temperature-independent vacuum energy, is divergent. 



However, this term cancels in (5.1), and thus f q (m) is finite. In figure |3]we plot the functions 
f q (m) for various q. These functions decay exponentially fast in m. 



Mm) 




m 



Figure 3: The Renyi entropy interpolating function f q (m) for a massive complex scalar 
field, plotted as a function of m for q = 1, 2, 3, 4, 5 (plots are darker for larger q). 
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5.2 Free massive Dirac fermions 



The calculation for free massive Dirac fermions proceeds analogously to that of the scalars. 



The action on H q is given in equation (4.1). The free energy J : q {m) can be found in equa- 



tion (4.15). Again, while the temperature-independent vacuum energy integral in equa- 



tion (4.15) is divergent, the function f q (m) is finite. In figure 4 we plot the Renyi entropy 



interpolating function f q {m) for various q. The function f q (m) decays exponentially fast in 



-0.05 



-0.10 - 



-0.15 



-0.20 




m 



Figure 4: The Renyi entropy interpolating function f q (m) for a massive Dirac field, plotted 
as a function of m for q = 1, 2, 3, 4, 5 (plots are darker for larger q). 



m and is stationary at m = for all q. 



6 Discussion 



In this paper we have presented detailed results for the Renyi entropies of d = 2 free conformal 
fields across a circle, using conformal mapping methods. Let us note that the approach using 



the conformal mapping to S 1 x M d may be easily generalized to d > 2 19,49 where the 
entanglement is measured across S . All we need in this approach is the density of states 
m H d , which may be found, for example, in |54|J 



In d = 3 this approach was used in [19 ,56 to calculate the Renyi entropies for a free 



conformal real scalar field (1.3). It is not hard to rederive this result using the formulae in 



our paper. As in section 3.2 the free energy T q is given by equation (3.22) with m = and 



For direct numerical calculations of the Renyi entropies in d > 3 see 55 . 
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a density of states T>{\) appropriate for H 3 (see for example 52 -54]): 



(6.1) 



Note that the density of states on M 3 is identical to that in flat space. The regulated volume 



of H 3 is calculated in 19 to be 



Voi(e 3 



-27rlog 



R 



(6.2) 



After regularization, we then find that the free energy JF q is simply given by 



T 



1 (R 
log 



360g 



(6.3) 



Upon substitution into (1.9) this leads to the Renyi entropies (1.3). 



This approach is also easily applied to a massless Weyl fermion in d — 3. The free 



energy T q of the theory on S 1 x H 3 is given by equation (4.15) with m = and density of 



states 53,54 



V{\) 



Voi(e 3 ) 



7T" 



X 2 + 



(6.4) 



This leads to the free energy 



~ 7 + 30g 2 , (R 
^ = l440^ 1Og 



(6.5) 



Note that this result agrees with a similar calculation in de Sitter space 57 . Substituting 

l + g)(7 + 37g 2 ) 



into (1.9), we find the Renyi entropy 

~S n = 



1440g 3 



log 



R 



(6.6) 



in agreement with 56 . The presence of the factor 1 + q in d = 3 is guaranteed by the fact 
that T q is expanded in odd powers of q. 

As a check of these results, we note that the universal part of the d = 3 entanglement 

7T 



entropy satisfies 58 



Si = — alo£ 



(6.7) 
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where a is the anomaly coefficient normalized in such a way that a = 1/90 for a real scalar, 
and a = 11/180 for a Weyl fermion. This general formula follows from the fact that the 
q — 1 calculation may be mapped to the partition function on S 4 whose logarithmic piece is 
determined by the Weyl anomaly. Another consistency check is that in the high temperature 
limit q — > 0, the curvature of H 3 may be neglected and S q must be proportional to the 
standard thermal entropy calculated on S 1 x IR 3 . As q — > 0, we find that S q /S q — > 7/4 which 
is indeed the correct ratio of thermal entropy of a Weyl fermion and a real scalar. 

Let us note that both the scalar and the fermion Renyi entropies may be expressed as 



where B = ^therm/^therm an d ^ = 180 fl — B, which follows from a free energy on S 1 x H 3 of 
the form 



B + (A - B)q 2 , _ (R 
360g 



" log T • (6-9) 



It is clear [49 1 that in d — 3 the Renyi entropies for g / 1 are not determined solely 
by the anomaly coefficient a, which is a protected quantity that is not corrected beyond 
one loop. We have seen that they also contain input from the thermal free energy, which 
in general receives various perturbative and non-perturbative corrections. Therefore, we do 
not expect S q for q ^ 1 to satisfy a "c-theorem." In particular, we do not expect them to 
be constant along lines of fixed points. For example, in the M = 4 SYM theory the Renyi 
entropies should depend on the 't Hooft coupling because the thermal free energy has such 



dependence 59 . The entanglement entropy Si, which is determined by a through (6.7) 



appears to be the only special quantity which is monotonic along RG flows and is stationary 



at the fixed points 60 61 . We believe that these statements are valid not only for d — 3, 
but for all dimensions d > 1. We hope to consider some d > 1 Renyi entropies in more detail 
in a later publication. 
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A Degeneracies on S* 3 and its g-fold branched coverings 
A.l Setup 



The parameterization (2.11) that yields (2.6) is closely related to the Hopf fibration from 



the case q — 1. The Hopf fibration is defined as 

Zx = cos <t e iM)/2 Z2 = sin l e i(-n^)/2 (A.l) 
2 2 

So after writing 9 = 9/2, r = (0 + ip)/2, and = (— + ip)/2 the metric becomes 

1 



ds 2 - - 
4 



d9 2 + sin 2 9 d(p 2 + (chj> + cos 9 dfj 



(A.2) 



which exhibits S 3 as an S 1 fibration over S 2 . The ranges of these angles are customarily 
taken to be < 9 < it, < < 2n, and < up < Att. 

To work with spinors it is convenient to consider the forms cij defined as 

crj + ia 2 = e i4,+iT [id9 + sin 9 cos 9(dr - d<p)} , a 3 = cos 2 9 dr + sin 2 9d(p . (A.3) 

These <7j are the left-invariant one forms on S 3 , and away from 9 = 0,tc/2 they are well- 
defined and non- vanishing. They satisfy 

d(Ji = e ijk aj A a k , (A A) 

and, in addition, the metric on S* 3 can be written as 

ds 2 = o\ + o\ + o\ . (A.5) 
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We will henceforth use the left-invariant forms dreibein. 

On the g-fold branched covering of S 3 (2.6) with q > 1 one can still use the forms a. 



given in (A. 3) as a dreibein. The only difference is that now the range of r is [0,27rg), with 
q taken to be a positive integer. 



A. 2 Degeneracies for fermions 

To study fermions, let us take the gamma matrices in 3d to be equal to the Pauli matrices: 



7i 



1 

1 0. 



72 



-i 

1 



73 



1 
-1 



(A.6) 



That the spinor ip is an eigenspinor of the Dirac operator with eigenvalue A means that 

1 



<7 a (^ + ^<7a^) =AV>- 



For solving (A. 7) we make the ansatz 



i/j = e 



im,T+in<f> 



>(0+t)/2 



/+(*) 
/-(*) 



The solution that is square-integrable at 9 = tt/2 can be written as 

= \ [(2m - 1) tantf - (2n - 1) cot 9} f{9) + f{9) 



where 



(cos9) m 2 (sin i 



2 2 



( m + n -A--)/(0) 



-Fi ( 5 — 2-, 5 — m + |; cos^ 6> 



(cos#) 



2 m (sini 



-m— n— A+4 



-m-n+A+| 



-m + §; cos 2 



(A.7) 



(Ai 



(A.9) 



if m > 
if m < 

(A.10) 



Imposing regularity at 9 = as well restricts the allowed values of A in terms of n and m. 
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In particular, there are several infinite series of eigenmodes: two with positive A, 



X = m + n+- + 2a, 

A = — m + n-\ h 2a . 

2 



and two with negative A: 



\ = -\ m + n+- + 2a l . 



A = — ( -m + n + - + 2a , . 



m > , 
m < 



n > a . 

~ 2 

1 

n > a . 



a G N 
a G N. 



m > , n > a, 

- ' - 2 



m<0, n > — — — a , 

A, 



a e N. 



a e N. 



(A.ll) 



(A.12) 



To completely solve the eigenvalue problem (A. 7) on the g-fold branched covering of S s 



we need to decide on the allowed values of n and m by requiring the spinor (A. 8) to be 



well-defined. Under either r — > T + 2irq or — > <p + 2tt we must have ip — >■ — ^, which implies 
e 2™ng _ Consequently, n G Z + | and m G - (Z + |). It can be proven by 



n 2ivin 



induction that the eigenvalues of the Dirac operator then are those given in eq. (4.3) with 
multiplicity (k + l)(k + 2) for each choice of sign. 



A.3 Vector fields 

We can also find a complete basis of vector fields, or equivalently one-forms, on S 3 . To do 
so, we find the eigenfunctions of the vector Laplacian 



A = *d*d + d*d* 



(A.13) 



and these eigenfunctions will be our basis. Hodge decomposition states that any one-form 
A can be written as 



A = B + dd> + H . 



d*B = 



AH = 



(A. 14) 



for some 0, B, and H. Since H 1 (S 3 ) = 0, there are no non-zero harmonic one-forms H on 
S 3 . So in order to find a basis of forms A, it is enough to find a basis of exact one-forms dip 
and a basis of co-closed one-forms B. 

The basis of exact one-forms is given by derivatives of the spherical harmonics. Indeed, 
if Y is a spherical harmonic satisfying V 2 Y = — *d*dY = —n(n + 2)Y, then dY is of course 
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exact and is an eigenf unction of A: 



AdY = d*d*dY = n(n + 2)dY . 



(A.15) 



Diagonalizing A in the subspace of co-closed one-forms is easier on a three-dimensional 
manifold than on a manifold of arbitrary dimension because only in three dimensions does 
the operator *d map one-forms to one- forms. In fact *d is a self-adjoint operator on i7 1 (5' 3 ). 
Suppose we find a co-closed B that is an eigenfunction of *d with eigenvalue A: 



d*B = , *dB = XB } 



(A.16) 



Then B is also an eigenfunction of A: 



AB = (*d) 2 B = X 2 B. 



(A.17) 



In order to solve the equations (A.16) we make the ansatz 



B = e imT+in<t> [e- i{T+ ^b + (6)o + + b (6)a 3 + e iir+ ^b_(9)a_] . (A.18) 
For any given A, the solution that is regular at 9 = tt /2 is: 



b±(6) = (m + riT A) 
b (9) = (m + n + X)(m + n- X)f(9) 



(A.19) 



where 



(1 , ,„ . / \m\ + n + X 2 + |m| +n - A 
/(6») = (cos6') |m| (sm6') n 2-Fi ( J — — , — 7; ; 1 + \m\ ;cos z 9 



(A.20) 



Requiring that (A.19) is regular at 9 = yields four towers of solutions: two with positive 
A: 



X = m — n + 2 + 2a , m > , n < a + 1 



A = 2 — m + n + 2a. 



m < 



n > a — 1 



(A.21) 



where a G N, and two similar towers with negative A. 

Counting carefully we find that the eigenvalues of *d in the space of co-closed one-forms 
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are of the form ±n, with degeneracies 

fn 2 -l ifneZ, n>2, 
g n ={ ~ (A.22) 

^m(m + l) if n = m + ^ , m G Z , m>l, 

for every integer p that satisfies 1 < p < q. 

B Useful mathematical formulae 

In this section we present some useful mathematical formulae. We begin with zeta function 
identities. For < a < 1 we have the identity 



C(z,a) = 



zir cos 2iran zir sin 2nan 



sin 



Euus Zi /i an a i\ \ — * sin z, , 
h cos — > - 

n=l n=l 



J2 log (\ + a 2 ] = 2 log [2 sinh(7rg |a|)] , 
h &(jf +a2 ^ = 21og[2cosh(vrg|a|)] . 



nez+2 



These sums follow from the more general formula 



(B.1) 



(27T) 1 - 

Taking derivatives at z — 0, — 1, — 2 gives 

; 1 cos 2n an 1 (2 log(27m) + 27 — 3) sin 2nan 

C (-2, a) = -— 2^ ~ 3 , 

n=l n=l 

. 1 *^sin27ron 1 (log(27rn) + 7 — 1) cos27ran 

C'(-l,a) = — V ^--^7^^ " L ^— L , B.2 

47r ^ n 2 2tc 2 ^ n 2 v ' 

n=l n=l 

. 1 cos27ran 1 (log(27m) +7)sin27ran 
C(0,a) = +-Z^ " " 

n=l n=l 

Two other useful identities are the regularized sums 



(B.3) 



^ log + + a 2^ = log [ 2 cosh(27rg |a|) - 2 cos(27ra)] . (B.4) 
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This relation in turn follows from the Poisson summation formula 

n + a 



— E/ 



1 



E' 



i2nka 



f(2nqk) 



(B.5) 



applied to 



/(w) = log (u; 2 + a 2 ) . 



(B.6) 



For t ^ one can simply calculate the inverse Fourier transform of /: 



/(*) 



dcu 
2^ 



e- lU]t log (u 2 + a 2 ) 



-\a\W 



t 



(B.7) 



The case t = requires special care because the expression for /(0) is divergent and requires 
regularization: 



/(0) 



do; 
2^ 



log (w 2 + a 2 ) 



d f duj 1 



a . 



(B.8) 



Using (B.6)-(B.8) one can show that eq. (B.5) reduces to eq. (B.4) 



C Chern-Simons theory on C q 

In this section we evaluate the free energy J^s of U{1) Chern-Simons (CS) theory on C q 
and reproduce the known result that J-cs is independent of q. The partition function of CS 
theory is 



^cs 



where Scs is the Chern-Simons action 

ik 



(c.i; 



S cs = — I AAdA 

47T 



(C.2) 



and ^ghost is an action involving ghosts required for proper gauge fixing. We use the gauge 
d * A = 0, in which case the ghost action is 



S g host = J d 3 xy^ [cV 2 c + 6VM J 



(C.3) 
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where c is a fermionic ghost and b is a bosonic ghost. The factor of 1/ ^/Vol(C 9 ) in the 
partition function is important in order for the answer to be independent of the metric on 
C q . The factor of 1/ (2tt) comes from the volume of the gauge group, which one customarily 
divides by. 



The partition function (C.I) yields the free energy 



JPfS = J-vec + ^-ghost + lQg(27r) + 1 log(g Vol(5 3 )) , (C.4) 

where we used the fact that Vol(C 9 ) = gVol^), and denoted by J-J ec the contribution from 
the Chern-Simons action ( |C2 ) in the gauge d*A = 0, and by J^& host the contribution from the 



ghosts. It is straightforward to show that F& host equals minus the free energy of a massless 
complex scalar on C g , while J 7 ™ can be computed from the eigenvalues and degeneracies in 



(A. 22). When q = 1 we have 

n=2 

^.ghost = I^( n + i)2 i ogn{n + 2 ) = 1 logvr + C(3) 



(C.5) 



4tt ' 

n=l 



which implies = | log fe. For the gth branched covering of S 3 we have: 



^ cs = ^ 1 cs + 1 ^ + EE-(-+ 1 ) 



<J— 1 OG 



p=l m=l 

9-1 



1 / p\ / p 



log ^ 2 — - - - log - 1 + - J ( /// I 



p=i 



(C.6) 

The fact that = J-^ s for all g can be understood from the fact that the topology of 



the branched covering C„ is the same as that of S 3 . The Chern-Simons partition function 



is a topological invariant, so it should indeed be independent of q. From eq. (1.9) it follows 
that all the Renyi entropies are equal: 



S<* = S° a = -\)o g k. (C.7) 
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